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Abstract 

A unified approach to energy-efficient power control is proposed for code-division multiple access 
(CDMA) networks. The approach is applicable to a large family of multiuser receivers including the 
matched filter, the decorrelator, the linear minimum mean-square error (MMSE) receiver, and the (nonlin- 
ear) optimal detectors. It exploits the linear relationship that has been shown to exist between the transmit 
power and the output signal-to-interference-plus-noise ratio (SIR) in the large-system limit. It is shown 
that, for this family of receivers, when users seek to selfishly maximize their own energy efficiency, the 
Nash equilibrium is SIR-balanced. In addition, a unified power control (UPC) algorithm for reaching the 
Nash equilibrium is proposed. The algorithm adjusts the user's transmit powers by iteratively computing 
the large-system multiuser efficiency, which is independent of instantaneous spreading sequences. The 
convergence of the algorithm is proved for the matched filter, the decorrelator, and the MMSE receiver, 
and is demonstrated by means of simulation for an optimal detector. Moreover, the performance of 
the algorithm in finite-size systems is studied and compared with that of a conventional power control 
scheme, in which user powers depend on the instantaneous spreading sequences. 
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Index Terms 

Code-division multiple access (CDMA), energy efficiency, game theory, large systems, multiuser 
detection, multiuser efficiency, Nash equilibrium, power control. 

I. Introduction 

Power control is used for interference management and resource allocation in both the down- 
link and the uplink of code-division multiple access (CDMA) networks [l]-[8]. In particular, 
in the uplink each user transmits just enough power to achieve the required quality of service 
(QoS) without causing excessive interference in the network. In recent years, power control for 
multiuser receivers has attracted much attention due to their superior performance as compared 
to the single-user matched filter. For example, power control algorithms for the linear minimum 
mean- square error (MMSE) receiver and successive interference cancellation receivers have been 
proposed in [4] and [7], respectively. In the proposed schemes, the output signal-to-interference- 
plus-noise ratio (SIR) of each is measured and then the user's transmit power is adjusted to 
achieve the desired SIR. 

Power control is conventionally modeled as an optimization problem under some quality 
of service (QoS) constraints. A practically appealing scheme is to minimize the total transmit 
power under the constraint that the output SIR of each user is above some lower bound. Under 
reasonable assumptions it is shown that the total transmit power is minimized when the SIR 
requirements are met with equality [1]. Another approach is to choose the transmit powers in 
such a way as to maximize the spectral efficiency (in bits/s/Hz) of the network, where the optimal 
strategy is essentially a water-filling scheme [6]. 

In recent years, game theory has been used to study power control in CDMA networks. In 
particular, power control can be modeled as a non-cooperative game in which each user selfishly 
maximizes its own utility. The strategy chosen by one user affects the performance of other users 
in the network through multiple-access interference (e.g., [9]— [15]). 

Reference [14] studies the cross-layer design problem of joint power control and multiuser 
detection using game theory. The utility function in this case measures the number of bits 
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transmitted per joule of energy consumed, which is particularly suitable for energy-constrained 
networks. Focusing on linear receivers, the Nash equilibrium for the non-cooperative game has 
been derived. A Nash equilibrium refers to a set of strategies for which no user can unilaterally 
improve its own utility [16]. It is also shown in [14] that for all linear receivers, the transmit 
powers of the users are SIR-balanced at Nash equilibrium (i.e., all users have the same output 
SIR). 

This work extends the results in [14] to a much larger family of receivers (including nonlinear 
ones) in the so-called large-system regime, where the number of users and the spreading factor 
are large with a given ratio. This is due to results in [17] where a linear relationship between 
the input power and the output SIR is shown to exist for a family of multiuser detectors in the 
large-system limit. Members of this family include well-known receivers such as the matched 
filter (MF), the decorrelator (DEC), the linear MMSE receiver, as well as the individually optimal 
(10) and jointly optimal (JO) multiuser detectors!] By exploiting the linear relationship, which is 
characterized by the multiuser efficiency, we propose a unified power control (UPC) algorithm 
for reaching the Nash equilibrium. The convergence of the proposed algorithm is proved for 
linear receivers and is demonstrated by means of simulation for an optimal detector. 

Since the UPC algorithm is based on large-system results, it does not depend on the instan- 
taneous spreading sequences. The true SIR, however, does depend on the spreading sequences. 
Therefore, the SIR achieved by the UPC algorithm deviates from the target SIR as the spreading 
sequence changes from one symbol to the next. This in turn will result in a loss in the utility. 
The performance of the UPC algorithm in finite-size systems is studied and it is shown that if 
the spreading factor is reasonably large, the SIR achieved by the UPC algorithm stays close to 
the target SIR most of the time which means that the loss is insignificant. 

The rest of the paper is organized as follows. Section [TT] provides the system model and relevant 
results in multiuser detection. Section [TIT] discusses the game-theoretic approach to power control. 

'The individually optimal detector minimizes the error probability of detecting an individual symbol whereas the jointly 
optimal detector minimizes the probability of error for detecting the entire vector of all users' symbols at a symbol interval 
[18]. The jointly optimal detector is often referred to as the maximum likelihood (ML) detector. 
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Fig. 1. (a) The DS-CDMA system with multiuser detection, where only the output for user k is shown, (b) The equivalent 
single-user Gaussian channel for user k; also shown is the companion channel. 

The UPC algorithm for reaching Nash equilibrium is proposed and studied in Section [IV] The 
performance of the UPC algorithm in finite-size systems is studied in Section |V] Simulation 
results are presented in Section |VI] before conclusions are drawn in Section IVIIl 



II. Multiuser Detection and Power Control 

A. System Model and Multiuser Detection 

Consider the uplink of a synchronous DS-CDMA system with K users and spreading factor 
N. Let p k , hk, and a 2 represent the transmit power, channel gain and the background noise 
variance (including other-cell interference) respectively, for user k. The received signal-to-noise 
ratio (SNR) for user k is then 

r 4 = ^. (i) 



The received signal (after chip-matched filtering) sampled at the chip rate over one symbol 
duration can be represented as 



A 



Y = sx + w = ^ v / rv^ + w, 



(2) 



k=l 



where Sfc and Xk are the spreading sequence and transmitted symbol of user k, respectively, 
and W ~ jV(0, I) is the noise vector consisting of independent standard Gaussian entries (see 
Fig. [Ha)). Random spreading sequence is assumed for all users, and the input symbols X k are 
assumed to be independent and identically distributed (i.i.d.) with unit variance. 



Suppose the receiver is a linear filter c k , the detection output can be represented as a sum of 
three independent components: 

X k = c T k Y (3) 
%X k + MAI, + V k (4) 



i.e., the desired signal, the MAI and Gaussian background noise. The quality of the output is 
effectively measured in terms of the output SIR, 

because the multiple-access interference (MAI fc ) is asymptotically Gaussian (e.g., [19]). The 
degradation in SNR due to the MAI is known as the multiuser efficiency, denoted by rj k . In 
fact, © can be regarded as an equivalent single-user Gaussian channel for user k as depicted 
in Fig. (2b), in which the function g(-) is an identity mapping, e.g., g(z) = z, Vz. 

B. Power Control 

We model power control as a non-cooperative game in which each user tries to maximize its 
own utility (see e.g., [9]— [15]). We follow [9] to define the utility of user k as 

T 

u k = — bits/joule, (6) 

Vk 

where T k is the throughput of user k, i.e., the net number of information bits delivered correctly 
per unit time (sometimes referred to as goodput). This utility function captures the tradeoff 
between throughput and battery life and is particularly suitable for applications where saving 
power is more important than achieving a high throughput. 
The throughput for user k can be quantified as 

T k = ^Rf(lk), (7) 

where L and M are the number of information bits and the total number of bits in a packet, 
respectively (i.e., L — M bits of overhead); R is the transmission rate and f(^ k ) is the efficiency 
function representing the packet success rate. The underlying assumption is that if a packet is 



in error, it will be retransmitted. The efficiency function, f(jk), is assumed to be increasing and 
S-shaped^ (sigmoidal) with /(oo) = 1. We also require that /(0) = to ensure that u k = when 
Pk = 0. These assumptions are valid in many practical systems (see [14] for further details). 
Combining © and ©, user k's utility is given by 

u k = — tj • (8) 

Mpk 

Using a sigmoidal efficiency function, it can be shown that the utility function in ([8]) is quasi- 
concave[] 

Suppose each user is allowed to vary its transmit power only in order to selfishly maximize 
the utility, the power control game is described as 

maxttfc for k = 1, • • ■ , K. (9) 

Pk 

For this non-cooperative game, a Nash equilibrium is a set of transmit powers (pi, . . . ,Pk) for 
which no user can unilaterally improve its own utility (i.e., a stable state). It has been shown in 
[14] that the following proposition holds for all linear receivers. 

Proposition 1: The utility-maximizing strategy for user k is given by p* k = min(pfc, P ma x), 
where pk is the transmit power that results in an output SIR equal to 7*, which is the solution 
to /(t) = 7/' (7)' an d Pmax is the maximum allowed power. Furthermore, the proposed power 
control game has a unique Nash equilibrium. 

This proposition implies that, at Nash equilibrium, the users' transmit powers are SIR-balanced. 
The key point in proving Proposition \T\ is that for all linear receivers there is a linear relationship 
between the output SIR and transmit power of a user. If we take the derivative of the utility 
function in © with respect to the transmit power and equate it to zero, we have 

Pk^filk) ~ f{lk) = 0. (10) 

OPk 

2 An increasing function is S-shaped if there is a point above which the function is strictly concave, and below which the 
function is strictly convex. 

3 A function is quasiconcave if there exists a point below which the function is non-decreasing, and above which the function 
is non-increasing. 



Based on the expression of the output SIR ©, we can write 

dp k Pk 

Therefore, the utility of user k is maximized when j k = 7*> the (positive) solution to 



(ID 



/(7) = lf{l) ■ (12) 

The existence of a Nash equilibrium is guaranteed because of the quasiconcavity of the utility 
function and its uniqueness is because of the one-to-one relationship between the user's transmit 
power and output SIR. 

III. Large Multiuser Systems and Power Control 

This paper studies power control in large systems. Mathematically, we consider the so-called 
large-system limit, where both the number of users and the spreading factor tend to infinity 
but with their ratio converging to a positive number, i.e., K/N — > a. Since, in a large system, 
the dependence between the SNRs Y k is reasonably weak, we assumed that Y k are i.i.d. with 
distribution P r at a given time. Moreover, P r varies slowly over time as the timescale of power 
control is much larger than the symbol interval. 

In general, the multiuser efficiency depends on the received SNRs, the spreading sequences as 
well as the type of detector. However, in the asymptotic case of large systems, the dependence on 
the spreading sequences vanishes and the received SNRs affect r] only through their distribution. 
In particular, the multiuser efficiency of the matched filter and the decorrelator are obtained as 

' 1 + aE{T} K } 

r) dec = l- a for a < 1 (14) 

while the efficiency of the (linear) MMSE receiver is the unique solution to the following fixed- 
point equation 

1 1 + oE\- -1 (15) 



where the expectation is over P r . 
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The analysis of the SIR using © does not directly apply to nonlinear receivers because the 
output of such receivers cannot be decomposed as a sum of the desired signal and independent 
interference; neither is the output asymptotically Gaussian. Remarkably, reference [17] finds that, 
under mild assumptions, the output of a nonlinear receiver converges in the large-system limit 
to a simple monotone function of a "hidden" Gaussian statistic conditioned on the input, i.e., 

X k ^g(y/V~ k X k + U k ) (16) 

where U k ~ A/"(0, rf~ 1 } is independent of X k . By applying an inverse of this function to 
the detection output X k , an equivalent conditionally Gaussian statistic Z k = ^/TlXk + U k 
is recovered. Each symbol X k traverses an equivalent single-user Gaussian channel, so that the 
output SIR (defined for the equivalent Gaussian statistic Z k ) completely characterizes the system 
performance. This result is referred to as the "decoupling principle." The equivalent channel is il- 
lustrated in Fig.QJb). Indeed, as far as the posterior probability P Xk \x k ls concerned, the multiuser 
model (Fig. [fla)) and the single-user model (Fig. QIb)) are asymptotically indistinguishable. 

A. Posterior Mean Estimators 

The decoupling principle is shown in [17] to hold for a broad family of multiuser receivers, 
called the posterior mean estimators (PME). Given the observation Y and the spreading matrix 
S, a PME computes the mean value of some posterior probability distribution gx|Y,s> which is 
conveniently denoted as 

(X) g = E 9 {X| Y,S} (17) 

where E g {-} stands for the expectation with respect to the measure q. 

In this work, the posterior gx|Y,s supplied to the PME is induced from the following postulated 
CDMA system, 

Y = SX' + gW (18) 

which differs from the actual channel © by only the input and the noise variance. In particular, 
the components of X' are i.i.d. with distribution qx, and the postulated noise level g serves as 



a control parameter. The posterior gx|Y,s is determined by q x and gY|x,s according to Bayes' 
formula 

/ i G x gx(x)g Y |x,s(y|x,S) 

9x|y,s(x y, S) = T — — — — . (19) 

J gx(x)qv|x,s(y|x, S) dx 

Indeed, the PME so defined is parameterized by (qx, g) and can be regarded as the optimal 
detector for the postulated multiuser system (PT8T) . In case the postulated posterior qx|Y,s is 
identical to px\Y,s, the PME is a soft version of the individually optimal detector. The postulated 
posterior, however, can also be chosen such that the PME becomes one of many other detectors, 
including but not limited to the matched filter, decorrelator, linear MMSE receiver, as well as 
the optimal detectors. Thus the concept of PME is generic and versatile. 

B. Decoupling Principle for PME 

Let qz\x,r-£ represent the input-output relationship of a scalar Gaussian channel, 

(20) 



-f (z-y/Tx^ 



qz\x,T;t(A x i V iC) = y^ ex P 2 
Similar to that in the multiuser setting, by postulating the input distribution to be qx, a posterior 
probability distribution qx\z,v& is induced from q x and qz\x,v-£ using Bayes' formula (cf. (fl9l). 
Thus we have a single-user companion channel defined by qx\z,T& which outputs a random 
variable X' given the channel output Z (Fig. djb)). A (generalized) single-user PME is defined 
naturally as: 

(X) q = E q {X\ Z,Y-i) (21) 

The probability law of the composite system depicted by Fig. [TJb) is determined by T, 77 and £. 

Proposition 2 ( [17]): Fix (a, Pr,Px, Qx, f?). The joint probability distribution of (X k , {X k ) q ) 
converges to the joint probability distribution of (X, g(\/T^X + U)) where U ~ A/"(0, r/ _1 ) and 

g(z)=E q {X\ Z = z,T;0, (22) 
where the multiuser efficiency r\ satisfies together with £ the coupled equations: 

r]- 1 = 1 +aK{T-£(T;r ] ,0}, (23a) 
r 1 = f? 2 + «E{r-V(r ; r / ,0}, (23b) 
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where the expectations are taken over P r . Here we define the mean squared error of the PME 



as 

, 2 



e{T;7i,Z)=El(x-(X) 
and also define the variance of the companion channel as 



I>,0, (24) 



V(T;ti,£)=e{(x'-(X) 



2 

g 



r;77,0. (25) 



In case of multiple solutions to (|23l ), (r],£) is chosen to minimize the free energy (see [17]). 

Proposition |2] reveals that, from an individual user's viewpoint, the input-output relationship 
of the multiuser channel concatenated with the multiuser receiver is asymptotically identical to 
that of the scalar Gaussian channel with a (nonlinear) decision function. The key performance 
measure is the equivalent SIR r\ of the scalar channel, which is easy to compute by solving 
the fixed-point equation (|23l) . since the functions £ and V can be easily computed given the 
distribution of the transmitted symbols and the type of receiver; and the expectations are taken 
with respect to the received SNR distribution. 

By choosing appropriate parameters (qx,p), the PME can be made to represent the matched 
filter, the decorrelator, and the linear MMSE detector as well as the individually and jointly 
optimal detectors. The multiuser efficiencies of the linear receivers are given as (fT3l) - (fT5l) . 
The multiuser efficiency of the individually optimal detector is found to satisfy the fixed-point 
equation 

— = 1 + oE\t -T e -jL tanh (rf°V - z y/rj^r) dz \ . (26) 

V™ [ J-oc V27T V /J 

This result generalizes to multirate systems as well [20]. 

In principle, the multiuser efficiency may be different for different users with different types 
of receivers. As the uplink is the focus of this work, it is assumed that all users use the same 
type of receiver and hance are subject to the same multiuser efficiency. 
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C. Power Control 

The decoupling principle implies that in large systems there is a linear relationship between 
the transmit power and the output SIR of each user, i.e., 



where the multiuser efficiency rj depends on the SNR distribution P r rather than the individual 
IVs. This is mainly due to the fact that in a large system, as one user's transmit power varies, 
the interference seen by the user essentially stays the same as long as the overall distribution of 
the received powers remains the same. 

Even though (|27l) is a large-system result, it is a very good approximation for most finite-size 
systems of practical interest. The linear relationship (|27l) implies that, for the family of PME 
receivers, (ITT]) is satisfied in large systems. Therefore, the argument given above for the linear 
receivers extends to a larger family of receivers that includes some nonlinear receivers such as 
the individually and joint optimal multiuser detectors. 

This means that, in the asymptotic case of large systems, the Nash equilibrium for the power 
control game in ([8]) is SIR-balanced for all the receivers belonging to the family of PME for 
which (1271) is satisfied. In other words, the Nash equilibrium is reached when each user transmits 
at a power level that achieves an output SIR equal to 7*, the solution of 7(7) = 7^(7) . It is 
interesting to note that the Nash equilibrium SIR, 7*, is independent of the type of receiver and 
depends only on physical-layer parameters such as modulation, coding and packet size. 

Clearly, any deviation of a user's output SIR from 7* results in a loss in the user's utility, or 
equivalently, its energy efficiency. For example, fixing other users' transmit powers, if the output 
SIR of user k is equal to 7 (instead of 7*), the resulting utility is given by 



Ik = V^k 



(27) 




(28) 



where u* k is the maximum utility for user k (corresponding to 7^ = 7*). 
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IV. Unified Power Control Algorithm 

In this section, we propose a unified power control algorithm for reaching the Nash equilibrium 
of the power control game, which is applicable to the family of PME receivers. The algorithm 
iteratively adjusts the transmit powers in order to reach an output SIR equal to 7*. While we 
assume balanced SIRs here as a property of the Nash equilibrium, the UPC algorithm is also 
applicable to the case of unequal target SIRs. 

The UPC Algorithm carries out the following iteration: 

1) Let n = 0, start with initial powers Pi(0), • • • ,p#(0). 

2) Use (|23T) and the power profile to compute the multiuser efficiency, 77 (n). 

3) For user k, update the powers according to 

^ n+ v=^i£- < 29) 

4) n=n+l, stop if convergence; otherwise, go to Step 2. 

In Step 2, while finding an analytical expression for the multiuser efficiency is difficult for 
most multiuser detectors, it can be easily obtained from (|23T) using numerical methods. Note 
that (1231) need to be solved only once per iteration for each user. The uplink receiver (e.g., 
base station) can, for example, compute r] and feed it back to the user terminal. The above 
algorithm is applicable to a large family of receivers which includes many popular receivers 
such as the matched filter, the decorrelator, and the MMSE detector as well as the individually 
and jointly optimal multiuser detectors. Note, for example, that for a matched filter receiver, the 
UPC algorithm becomes the same as the bilinear power control algorithm proposed in [21] for 
minimizing the SIR error. In fact, one may also extend Proposition [2] to the case that different 
users use different type of receivers and Algorithm 1 will apply by replacing i] by individual 
efficiencies rjk. 

It is apparently a dilemma that the success of the power control scheme depends on the 
assumption that the SNR distribution is fixed, while the purpose of power control is to adjust 
the powers which may affect the SNR distribution. This, however, can be resolved naturally 
in practice, where the number of users is finite, by an iterative power control algorithm and 
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by replacing the expectations in (|23l) with an average over all users' received SNRs (or their 
estimates). For example, (fT51) can be expressed as 

i = 1 + fETT i F (30) 
V K ^ 1 + rjV k 

At any rate, the UPC scheme provides a viable mechanism to approach the Nash equilibrium. 

There are certainly other practical issues that need further investigation. For example, the UPC 

algorithm requires the background noise and other-cell interference power to be measured. This 

could be done by silencing the users in a cell in a coordinated manner. In addition, while we 

will show in Section [VI] that the UPC algorithm converges very quickly, the convergence rate 

and effects of estimation/measurement errors on the convergence of the algorithm require further 

analysis. 

We now prove the convergence of the UPC algorithm for the matched filter, the decorrelator, 
and the MMSE receiverQ Let T = [Tj, • • • , IV] and define an interference function, 

m = • < 3 » 

The dependence of the multiuser efficiency on T is explicitly shown here. By (|29l) and (1311 ), the 
UPC algorithm can be expressed as 

r fc (n + l) = /(r(n)), k = l,...,K. (32) 



Proposition 3: For the matched filter, the decorrelator, and the MMSE receiver, if there exists 
a f such that f\ > I(T), k = 1, . . . , K, then for every initial vector T(0), the recursion (|32l) 
converges to the unique fixed point solution of T* k = I(T*), k = 1, . . . , K. Furthermore, for any 
feasible f (i.e., f fe > /(f) for all k), T* < f fe for all k. 

Proof: The existence of a T implies that a feasible SNR vector exists for achieving 7*. It 
suffices then to show that I(T) is a standard interference function [3], i.e., for all T with > 
for all k, the following three properties are satisfied: 1) Positivity: I(T) > 0; 2) Monotonicity: If 

4 The convergence analysis for a general receiver remains an open problem because l !23t is difficult to work with. 
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T'fc > T k for all k, then I(T') > 1(F); 3) Scalability: For all 9 > 1, 0I(r) > I(6T). Evidently, 
it is equivalent to showing the three properties for I(T) = l/r](T). 
Positivity of I(T) is trivial by (I23al) for all receivers since 77 G [0, 1]. 

Consider first the matched filter. The multiuser efficiency is given by (fl"3T) . where E{T} = 

Ef=i r fe/^- If r ' > r > then E { r '} > E { r ) and hence ^( r ') > ^( r )- To P rove tne tnird 
property, note that, for 9 > 1, /(0r) = 1 + aE{9T} <9 + a9E{T} = 91(F). 

Consider now the decorrelator, the multiuser efficiency of which is constant 77 = 1 — a > 
for a < 1. Properties 2 and 3 are trivial. 

For the MMSE receiver, the multiuser efficiency is the solution to (TT3T ). or equivalently, the 
unique solution of 

77 + «E 1^^-1 = 1. (33) 

Note that the left-hand side of (|33l increases if both 77 and T increase. Thus if T' > T, we must 
have r](T') < r)(T) to maintain the equality. Hence, I(r') > I(T). To prove the third property, 
let us define 7/ = i](9T) and 77" = #77', where 9 > 1. Therefore, 

77" + a0E j^^-j =9. (34) 

Showing 91(F) > I(9F) is equivalent to showing 77" > 77. Since the left-hand side of (|33l) is 
increasing in 77, and 

ji , „ TT? J 1 i 1 1 11 1 l 



"" + « E {^TT}- 1 H 1 -5j" >x ' <35) 

we must have 77" > 77. Therefore, 9I(T) > I(9T). This completes the proof. 



V. Performance Evaluation and Discussion 



The UPC algorithm described in Section [TV] is a large-system approach and hence independent 
of the spreading sequences. Therefore, after convergence, the transmit powers need not be updated 
as long as the channel gains remain static. Evidently, the actual SIRs depend on the spreading 
sequences (see for example ©). Consequently, as the spreading sequence changes, the output 
SIRs achieved by the UPC algorithm fluctuate around the target SIR. This fluctuation results in 
a loss in energy efficiency (see (1281)). 
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The question of interest is: if the UPC algorithm is used, how close will the SIRs be to the 
target SIR? In the following, we study the performance of the UPC algorithm for finite-size 
systems and compare it with that of an SIR-based algorithm. We focus on the decorrelator and 
the linear MMSE receiver. 

A. The Decorrelator 

For the decorrelator, it is sensible to assume a < 1. The large-system multiuser efficiency is 
given by rf ec — 1 — a. Hence the SNR dictated by the UPC algorithm is 

Tl = r dec =-^—, k=l,---,K, (36) 
I — a 

which should lead to an SIR of 7* in the large-system limit. 

However, the actual output SIR for a finite-size system is given by 



Ik 



7 



1 — a 



kk 



(37) 



where S = [s 1; s 2 , • • • , s K ] and [ • ] kk extracts the kth diagonal entry of a matrix. 

It has been shown that in large systems, the distribution of 1/ [(S 2 ^)" 1 ] kk can be approximated 
by a beta distribution with parameters (N — K + 1, K — 1) [22]. As a result, the probability 
density function of 7^ is given approximately by 

7 VrSeJ 5(^-^ + 1,^-1) 1 j 

where z < T* dec and B(a, b) = J* t a ^ 1 (l — t) b ~ l dt. Therefore, as the spreading sequences change 
from symbol to symbol, the probability that 7^ stays within A dB of 7* is given by 

PAAec = Pr {|7de C (dB) - 7 *(dB)| < A} (39) 



i 



f, d Jz)dz , (40) 



A A 

where 7^ = 10~Tf>7* and = IO107*. 



Alternatively, the fluctuation of the actual SIR around 7* can be approximated less accurately 
by a Gaussian distribution with variance [23] 
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i.e., ^ dec ~ Af (7*, C 2 ) • Therefore, the probability that 7*. stays within A dB of 7* is approximately 
given by 

PlTc = $ f 2 ^-) " $ (^^P") ' (42) 

where $(■) is the cumulative distribution function of the standard Gaussian distribution. 

B. The MMSE Receiver 

If the linear MMSE receiver is used and all users have the same target SIR, 7*, the steady-state 
SNRs will be identical to T* = 7*/r/ after the UPC algorithm converges, where the multiuser 
efficiency is given by 



!-« 1 1 //-, ^2 2(1 + a) 1 

It can be shown that the fluctuation of the true SIR around 7* is approximately Gaussian with 
variance [23], [24]: 



C . = ^_i2L rTI . (44) 



1=2 

The probability that 7 fc stays within A dB of 7* admits a similar expression to (1421) using the 
function $(•). 

It is seen from these approximations that the variance of fluctuations of SIR decreases as 
1/N. In the following section, we demonstrate the performance of the UPC algorithm using 
simulations and also investigate the accuracy of the theoretical approximations discussed above. 

VI. Numerical Results 

Consider the uplink of a randomly spread DS-CDMA system with K users and spreading 
factor N. The noise variance is assumed to be a 2 = 1.6 x 1CT 14 . We use 7(7) = (1 — e~ 7 ) M as 
the efficiency functions with 7* = 6.4 (=8.1 dB). 

We first demonstrate the convergence of the UPC algorithm assuming K = 8 and N = 32. 
The channel gain for user k is given by hk = 0.1 x where dk is the distance of user k from 

5 This is a useful example for the efficiency function and serves as an approximation to the packet success rate that is very 
reasonable for moderate to large packet sizes. 
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Fig. 2. Transmit powers for the ML, MMSE, and decorrelator, using the UPC algorithm (N = 32 and K = 8). 

the uplink receiver (e.g., base station). Assume dk = 100 + 10A; in meters. We implement the 
UPC algorithm for the decorrelator and the MMSE receiver as well as the maximum likelihood 
detector. 

Fig. [2] plots the transmit powers for users 1, 4, and 8 at the end of each iteration. It is seen 
that for all three receiver types, the UPC algorithm converges quickly to steady- state values. 
The results are similar when the initial power values and/or K and N are changed. It is also 
observed that the steady-state transmit powers for the decorrelator and the MMSE receiver are 
close to those of the ML detector (in this case, the difference is less than 22%). This means that 
in terms of energy efficiency, which is quantified by the utility achieved at Nash equilibrium, 
the decorrelator and the MMSE receiver are almost as good as the ML detector. 

We next investigate the fluctuation of the SIR and bit-error-rate (BER) achieved by the 
(large-system) UPC algorithm against perfect power control where the SIR is computed using 
instantaneous spreading sequences (labeled SIR-based in plots). Fig. [3] shows the SIR and bit- 
error-rate (BER) of user 1 using the MMSE detector. It is seen that the SIR-based algorithm 
achieves the target SIR, 7*, at all time whereas the output SIR for the UPC algorithm fluctuates 



18 



9 Q 



V 



-O- UPC 
SIR-based 



15 20 25 30 35 40 

Iteration 



88 « , 



3 t> e ,9 



-O- UPC 
x SlR-ba 
9 



15 20 



30 35 



40 45 



Fig. 3. User 1 output SIR and BER for the UPC algorithm and SIR-based algorithm with the MMSE receiver (N = 32 and 
K = 8). 



around the target SIR as the spreading sequences change. Also, the fluctuations in the BER are 
larger when the UPC algorithm is used. 

To evaluate the accuracy of the theoretical approximations given in Section |Vl Figs. |4] and [5] 
plot the cumulative distribution function (CDF) of the output SIR 7 for the decorrelator and 
the MMSE receiver with different spreading factors and both low and high system loads. The 
plots show CDFs obtained from simulation (based on 100,000 realizations) as well as those 
predicted by the theoretical approximations given in Section |V] It is seen from the figures that 
the theoretical approximations become more accurate as the spreading factor increases. Also, 
the approximations are generally more accurate when the system load is low. Note that, for the 
decorrelator, the approximation based on a beta distribution is slightly more accurate than the 
one based on a Gaussian distribution, especially for small spreading factors and large system 
loads. 

To quantify the discrepancies between the simulation results and the theoretical approxima- 
tions, we compute P^dec an d Pa,mmse using the CDFs obtained from simulations as well as 
those predicted by theory (see (|39| ) and (l42l)). Table [I] shows the results for different spreading 
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Fig. 4. CDFs of 7 for the decorrelator. 




Fig. 5. CDFs of 7 for the MMSE receiver. 
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TABLE I 

Summary of results for the decorrelator and the MMSE receiver 



N 


rjsim 
" ldB, dec 

a = 0.25 


rjbeta 
MdB, dec 

a = 0.25 


Tjnorm 
MdB, dec 

a = 0.25 


rjsim 
-MdB, dec 

a = 0.75 


rjbeta 
MdB, dec 

a = 0.75 


Tjnorm 
MdB, dec 

a = 0.75 


Tjsim 

MdB, MMSE 
a = 0.25 


Tjnorm 
MdB, A/ MSE 

a = 0.25 


rjsim 

MdB, MMSE 

a = 0.75 


pnorm 
MdB, A/ MSE 

a = 0.75 


16 


0.77 


0.87 


0.74 


0.28 


0.19 


0.30 


0.93 


0.46 


0.41 


0.33 


64 


0.98 


1.0 


0.97 


0.54 


0.64 


0.55 


0.99 


0.76 


0.74 


0.61 


256 


1.0 


1.0 


1.0 


0.87 


0.96 


0.87 


1.0 


0.98 


0.98 


0.91 



factors and system loads for A = 1 dB. The numbers in the table represent the probability that 7 
is within 1 dB of 7*. Based on (l28~l) . a 1-dB increase in the output SIR results in 10% loss in the 
user's utility. The probabilities obtained by simulation suggest that the UPC algorithm performs 
better for the MMSE receiver than for the decorrelator. It is also seen from the table that when 
the spreading factor is small, the fluctuations in the output SIR are considerable, especially 
when the system load in high. The performance improves as the spreading factor increases. For 
example, for the MMSE receiver, when N = 256 and a = 0.75, the SIR stays within 1 dB of 
the target SIR 98% of the time. It is also observed that the theoretical approximations for the 
MMSE detector are pessimistic. For the decorrelator, while approximating the SIR by a beta 
distribution is more accurate (see Fig. HI), the values obtained for Pa,mmse by the Gaussian 
approximation are closer to the simulation results. This is because the slope of the CDF of 7 
is closer to the slope of the Gaussian CDF. Since Pa,mmse heavily depends on the slope of 
the CDF, it is more accurately predicted by the Gaussian approximation (rather than the beta 
approximation). 

VII. Conclusion 

A unified approach to energy-efficient power control in large systems is proposed, which is 
applicable to a large family of linear and nonlinear multiuser receivers. The approach exploits 
the linear relationship between the transmit power and the output SIR in large systems. Taking 
a non-cooperative game-theoretic approach with emphasis on energy efficiency, it is shown 
that the Nash equilibrium is SIR-balanced not only for linear receivers but also for some 
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nonlinear receivers such as the individually and jointly optimal multiuser detectors. In addition, 
a unified power control algorithm for reaching the Nash equilibrium is proposed. It would be 
straightforward to extend the unified approach to multirate and multicarrier systems based on 
related large-system results [20], [25]. 
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